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We discuss the semiclassical limit of the entanglement for the class of closed pure systems. By 
means of analytical and numerical calculations we obtain two main results: (i) the short-time 
entanglement does not depend on Planck's constant and (ii) the long-time entanglement increases 
as more semiclassical regimes are attained. On one hand, this result is in contrast with the idea that 
the entanglement should be destroyed when the macroscopic limit is reached. On the other hand, 
it emphasizes the role played by decoherence in the process of emergence of the classical world. We 
also found that, for Gaussian initial states, the entanglement dynamics may be described by an 
entirely classical entropy in the semiclassical limit. 
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I. INTRODUCTION 

Understanding the emergence of all aspects of the clas- 
sical world from the quantum theory is far from being a 
trivial task. In fact, it is a recurrent question which has 
been asked since the early days of quantum mechanics. 
However, at the level of the expectation values of position 
and momentum, the classical limit is reasonably under- 
stood 0, 0, [3] , according to the conceptual scheme we 
try to sketch in the next few lines. 

Concerning the Newtonian trajectories, which are 
based on the assumption of complete knowledge of the 
initial conditions, their emergence from the quantum ex- 
pectation values is predicted by the Ehrcnfcst theorem, 
valid for initially localized states. Within the Ehrenfest 
time scale, expectation values agree with the classical 
trajectories and, in this sense, the quantum mechanics 
recovers the determinism of the Newtonian theory. The 
classical limit (h — ► 0) in this case is in general asymp- 
totic and does not require any further mechanism such as 
decoherence. In other words, there exists a formal clas- 
sical limit for expectation values even in the context of 
closed systems. Although such a mathematical limit is 
not unique, one possible way to obtain it is particularly 
convenient: a two-degree of freedom classical Hamilto- 
nian may be produced with the assumption of a separa- 
ble dynamics of coherent states @ . In this approach, the 
classical theory emerges as a consequence of the absence 
of the entanglement. 

After the Ehrenfest time scale, the Newtonian mechan- 
ics is no longer able to mimic the quantum expectation 
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values of observables and the single trajectory determin- 
ism is missed. It means that predictions can only be 
made by means of averages over ensembles. But, even 
in this case, classical equations of motion are applicable 
(within the Liouville formalism). This is ensured by a 
mechanism known as decoherence [(|. Environment, to 
which all systems are inevitably coupled, destroys quan- 
tum correlations and turns quantum predictions identi- 
cal to those offered by the classical statistical formalism. 
Since the decoherence time is in general much smaller 
than the Ehrenfest time, the correspondence between 
quantum and classical world will be complete for all in- 
stants of the dynamics. 

Within the above described lines of thought, a natural 
step before the introduction of the action of the environ- 
ment is to analyze the classical limit of an intrinsically 
quantum quantity - the entanglement - for the class of 
closed pure systems, for which the theoretical semiclassi- 
cal limit (h — ► 0) is well established at the level of expec- 
tation values. Here, a word concerning the entanglement 
is in order, recalling that the entanglement is the great- 
est resource offered "exclusively" by the quantum world, 
which allows for the realization of the most challenging 
ideas, like quantum information processing and quantum 
computation 0,0, HI- Since entanglement is regarded as 
an intrinsically quantum property with no classical ana- 
log, it seems to be an adequate modern tool for the study 
of the semiclassical behavior of quantum mechanics. 

We then ask about the behavior of the entanglement in 
such situations in which an underlying classical dynam- 
ics is well-known. This question motivates us to analyze 
the semiclassical limit of the entanglement dynamics for 
closed pure systems in which classicality is achieved by 
means of macroscopic coherent states, i.e., minimal un- 
certainty states | a) with \a\ 2 S> 1. 

At a first sight, one may be induced to expect a de- 
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crease in the entanglement in the semiclassical limit. 
However, we show analytically and also by means of nu- 
merical examples that this is not a precise idea: in fact, 
the short-time entanglement is /i-independent for closed 
systems (section II). 

This paper is organized as follows: in section II we dis- 
cuss a short time expansion for the entanglement in an 
arbitrary pure bipartite system. In Section III a semi- 
classical analysis of the whole dynamics of the entangle- 
ment is numerically investigated for two different mod- 
els and an explanation in terms of mechanisms that are 
classically understandable as well is given. Section IV is 
reserved for some concluding remarks. 



II. THE SHORT-TIME ENTANGLEMENT 

For globally pure bipartite systems there exist suitable 
cntropic quantities obeying all requirements for an ade- 
quate entanglement measure jlOt 1 1 1| . Particularly useful 
is the reduced linear entropy, which offers in general the 
same informations given by the von Neumann entropy 
but with a much less computational effort. 

The reduced linear entropy (RLE) for a bipartite sys- 
tem belonging to the Hilbcrt sps.cc & — & u S v is defined 
by 



S(t) = l-Tr z [p z 2 (t)}, 



(1) 



in which. p u ,v 

(t) = Tr v>u [p(t)} is the reduced density op- 
erator. The index z denotes the subsystems "u" and "v" . 
The density operator p satisfies the von Neumann equa- 
tion ihp = [H,p], being H the Hamiltonian operator. 
This is a dynamical measure of the purity of the sub- 
system z in the unitary bipartite dynamics of initially 
pure states. For this case, the Schmidt decomposition 
[l2j guarantees that RLE is also a measure of the en- 
tanglement between the parts |l3[ . Furthermore, the 
Araki-Lieb inequality 01 guarantees that the RLE of 
the subsystems are identical for pure states. 

In order to obtain a semiclassical expansion for the en- 
tanglement in the approximation of short times we write 
the Dyson's series for the density operator up to the order 
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where pa = \uo)(uo\ ® |wo)( w o| is an initially disentangled 
pure state composed by the product of bosonic coherent 
states corresponding to the subsystems u and v. Now 
we apply the partial trace over the subsystem u in the 
coherent state basis to get 



Pv(t) 



P i 0) 



1 



t + o Pv 



(2) ,2 



(3a) 



P v " } = \vo)(vo 

(i) / d 2 ui (in\[H, pojlm) 



pi 2) = 



d 2 u x ( Ul \[H, [H,p }}\ Ul ) 



(3b) 
(3c) 

(3d) 



TT {iK) 2 

Then, the RLE defined by JTJ may be written as 

S(t)~S^ + SWt + SWt 2 , (4a) 



S<® = l-Tr v 
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where Ti: [.] = /^<t;i|.|wi). 

Straightforward manipulations on this equations yield 
5(0) = 5(1) = g Consequently, S(t) ~ S^t 2 . This 
result, which has already appeared in literature ^5j, is a 
direct consequence of both the purity and the separability 
of the initial state. 

Defining the dimcnsionless Hamiltonian 



Ht 
H=— , 
n 

we put (|4a|) in the following compact form: 



(5) 



S(t)~2 Coo(H)+Cn(H)-Cio(H)-Coi(H) , (6) 



in which he have defined the correlations 

Coo = ((V'o|H|^o)) 2 , 

Coi = (^|Hpi(0)H|Vto>. 
Cio = (^ |Hp2(0)H|Vo), 
C n = (V-o|H 2 |^o), 



(7a) 

(7b) 
(7c) 
(7d) 



with \ipo) — l u o) ® bo)- Interestingly, by © we may 
put the conditions for the existence of short-time entan- 
glement in terms of the inequality 



2 C 



Cix — Ci 



Coi) > 0. 



(8) 



Now we look for a semiclassical expansion for the short- 
time entanglement given by 10 in the basis of coherent 
states. We start by defining a general two-degree of free- 
dom classical Hamiltonian. 

T-C(q u ,p u ,qv,Pv) = ^2 C nmlkqlPu 1vP k v ( 9 ) 

This is an explicit fi-independent function. Following the 
ordered quantization procedure we obtain the corre- 
sponding Hamiltonian operator 

H=S U S V CrunlkQZPuQlPv, (10) 

n.m.l.k 
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with [Q Z ,P Z ] — ih. S z = exp (— ^dQ z dp z ) is the sym- 
metric ordering operator £|. This is a suitable quan- 
tization method for our purposes since it factorizes the 
dependences in h. There are other sources of H, namely, 
the operator H and the parameterizations of the coherent 
states for the phase space. 

The calculations of all terms in © are made by in- 
serting HI U| l in 0. Then, using the unities of the co- 
herent states basis, namely, l u = J ^-^\ui)(ui\ and 
l v = j ^-p-\vi)(vi\, we are lead to 



2t 2 



S(t) 
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tt{u v \H\u v{){u q vi \H\u Vq) , (11) 



The function (v,qVq\H\uiVi) is the most general kernel we 
have to manipulate in the calculation of the correlations. 
By <|10[) we obtain 

(u V \H\uiVl) = ^ c nmkl(uo\SuQZP™\ui) 

x{w\S v Q k v P l v \ Vl ). (12) 

The properties of the symmetric ordering operator Q 
allow us to write 



( Uo \s u Q^p: n \ Ul ) i du „ du 

— ■ — = e 2 o 
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with similar expressions for v. Notice that q lJ and p lJ 
correspond to the usual canonical phase space pair when 
i = j. Otherwise, when i ^ j, they are complex numbers 
composed by combinations of distinct canonical pairs. 

Given the usual parametrization of the coherent state 
label u (or v) for the phase space (Eq.lJS}) we may write 

® u = v§ ~ l ^P" ) ' Considering smooth wave packets 
at t = 0, we may regard this partial derivative as a small 
parameter and the exponential operator in l|13a|) may be 
expanded in few terms. After performing such adequate 
expansions in (|13af) we return to l|12|) to calculate the 
semiclassical expansion of the kernel. 

The next step is to calculate the integrals in <| 1 1 ft .which 
is done by means of the following useful formula [J] : 



\(u \u}\ 2 f(u ,u) = 

TT 



U — Uq 



(14) 



with = + d 2 u . This relation emerges from the 
fact that |(mo|m)| 2 is a Gaussian weight function. Once 



again we may get a semiclassical expression by expanding 
the exponential operator in the first orders in h. In fact, 
the procedure described up to here must be performed in 
a consistent way such that the short-time RLE may be 
written as 



S(t) ~ t 2 



o(h- 2 ) + oih- 1 ) + o{h°) + 



(15) 



Notice that the classical limit h — > seems to yield a di- 
vergence of the short-time entanglement, what would be 
a quite unexpected result. However, an exhaustive calcu- 
lation shows that 0(/r 2 ) = O^' 1 ) = and O(h ) ^ 0, 
in general. This derivation is really lengthy and tedious, 
but does not require any further recipe beyond that one 
presented above and it will be omitted here. 

Concluding this calculation, in a semiclassical regime, 
where h is finite, but arbitrarily smaller than a typical 
action of the system, the RLE is given by 



s{t) ~ t 2 o{h°), 



(16) 



showing that the short-time entanglement does not de- 
pend on h in a semiclassical regime. 

The constant O(H ) is constituted by a sum of terms 
involving first-order derivatives (in phase space variables) 
of the Hamiltonian function 10. Note that, the possibil- 
ity of expressing the first contributions to the entangle- 
ment dynamics in terms of classical quantities is a direct 
consequence of the special basis used in the calculation. 

Result (|16f) . which is general within the class of bi- 
partite bosonic systems in pure states, attests that the 
entanglement does not tend asymptotically to the classi- 
cally expected limit S(t) — 0, as does those well-behaved 
quantities like expectation values of canonical operators 
Q and P. 

Also it has to be remarked that the quadratic depen- 
dence in time found for the short-time entanglement is 
independent of the integrability of the underlying classi- 
cal system, i.e., it is always algebraic no matter whether 
the classical system is chaotic or regular. In an ap- 
proach based on the Loschmidt echo dynamics |l6l | it has 
been predicted different behavior for the entanglement 
in chaotic and regular regimes. However, this has not 
been done for the time scale we are considering, and fur- 
thermore, in that case the entanglement generation is ob- 
tained by means of two distinct unitary evolutions (echo 
operator). In this sense, there is no conflict with our re- 
sults. In fact, the quadratic dependence of the short-time 
linear entropy has already been found in [l5j . It is just 
a consequence of the initial separability of the quantum 
state. Recent studies also confirm our results concern- 
ing the algebraic dependence [T^ | and the ^-independence 
[l8j of the short-time entanglement. 

Next we present two numerical examples confirming 
the above result and also discuss the long-time behavior 
of the entanglement . 
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III. NUMERICAL ANALYSIS 

A. Classical States 

Since wc arc interested in the quantum-to-classical as- 
pect of the entanglement, we avoid initial states that are 
initially entangled, mixed or delocalizcd. The natural 
choice is the minimum uncertainty coherent states de- 
fined by 



\v) = D(v) \r) 



cxp [va — v a\, 



D(v) 



cxp 



arctan \ v\ 



(vJ+ -u*J_) 



(17a) 



(17b) 



D was defined respectively for the harmonic oscillator 
and for the angular momentum (spins). The reference 
state \r) stands for the ground states |0), in the Fock 
basis, and \J,—J), in the angular momentum basis. 

The coherent state label, v, has the following common 
paramctrization in terms of classical phase space vari- 
ables: 



■in 



q + ip 



(18) 



y/hJ-(q*+p>) 



For bosonic coherent states, the classicality is attained 
by taking |w| 2 sufficiently macroscopic to guarantee that 
the average energy (v\H\v) be much greater than a typ- 
ical spectral distance E n+ i — E n . In this case, tak- 
ing \v\ 2 — > oo is mathematically equivalent to requir- 
ing h — > 0. On the other hand, for spin coherent states, 
the classicality emerges by means of a more sophisticated 
limit 19]: J ->• oo, K -> and h^J J{ J + 1) w hJ = 1. 

We consider Hamiltonian operators like 



H = E x <g> 1 2 + li ® H 2 + H 12 , 



(19) 



the last term being a nonlinear coupling which is capable 
to entangle initially separable coherent states. 



B. Dicke Model 

First, we investigate the entanglement in one of 
the most important models in quantum optics (QO) 
and more recently in quantum computation (QC), 
namely, the TV-atom Jaynes-Cummings Model (or Tavis- 
Cummings model) |2dl |. In QO it is mostly known in 
the so called rotating wave approximation (RWA), and 
describes the coupling of a monochromatic field of fre- 
quency oj with N non-interacting two-level atoms with 
level separation e. Recently in connection with QC, the 



model received renewed attention from solid-state com- 
munity working on 'phonon cavity quantum dynamics' 
pll [23, Josephson Junctions and quantum dots [2^, 
quantum chaos |24j and quantum phase transitions |25| . 
These works are interested on the model in the origi- 
nal form conceived by Dicke [26| without the RWA. The 
Hamiltonian in the last case is composed by the following 
terms 



Hi = e J z , H 2 = hioa/i 



H12 



G 



2. J 



(aJ+ + a f J_ 



G 
s/2J 



(a t J+ + aJ_) , 



(20) 



in which the operators a and a) are the usual harmonic 
oscillator creation and destruction operators associated 
to the field, and J+, J_, J z are angular momentum ladder 
operators. Here, the spin algebra is associated to the 
atoms, being the total spin given by J = N/2. G and G' 
are real coupling constants, and generally could be taken 
unequal. 

This model possesses a classical counterpart H N . 1C u(J) 
obtained in ref. |2]j which presents several interesting 
properties, e.g., chaotic behavior. Another peculiar fea- 
ture is the scaling property 



TL-, 



,(J) H 

JV-JCM {J') 



J 



J 1 



(21) 



If this relation is satisfied for the pairs [E(J),J] and 
[E(J'), J'], then the associated dynamics will be totally 
equivalent. 

The initial state = \v a ) <8 \t)f) is thus constructed 
by means of the following process. We use the formula 
rj = ri \fj in order to re-scale the vector of initial con- 
ditions rj = (q a ,Pa, Qf,Pf), keeping the dynamics for 
,/ = 1 as a reference. This guarantees automatically that 
we obtain a dynamics [E(J), J] equivalent to the unitary 
[E(l), 1]. Then, wc construct the coherent initial state by 
putting Vj in (|T%)) and (|T7Jl . Resuming, we chose initial 
quantum states with different J, but producing always 
the same classical dynamics. The scmiclassical condi- 
tions in this case are satisfied by increasing J and keep- 
ing h constant (as it indeed occurs in nature). It may be 
shown that this is totally equivalent to the mathemat- 
ical condition h = 1/J, with J — * 00 and fixed initial 
conditions. 

In Fig^ we show numerical results for the entangle- 
ment as a function of time for several values of J. Ac- 
tually, for long times, the quantity plotted is a kind of 
mean entanglement of the states 'associated' to a given 
classical trajectory. We calculated it as follows: given a 
certain classical trajectory, we choose a set of M initial 
states iV'i(O)) centered at points along the trajectory, and 
calculate the respective RLE's, Si(t). Then, the averaged 
quantity is given by 



S m (t) 



1 M 



(22) 
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Such a calculation allows us to observe a smooth mean 
behavior for entanglement without the characteristic os- 
cillations associated to the border effect j^. The mean 
long-time entanglement is then suitably described by a 
fitting expression [2^ given by 



S m (t) = A Q (l-e- Alt ), 
where A and A\ are fitting parameters. 

1.0 i ■ , ■ , ■ , — 



(23) 
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(Gaussian) states. Such a ^-independence induces us to 
believe that the rising in the entanglement is determined 
just by classical sources, as the local departure of neigh- 
boring phase space points. This is indeed suggested by 
previous studies on quantum chaos |3C( . But the main 
point is that entanglement does not diminish as more 
macroscopic regimes are attained. 

The second aspect may seem even more curious at first 
sight, since it means that the entanglement increases as 
we tend to a more semiclassical limit. We will re-take 
this analysis after our second example. 



C. Coupled Nonlinear Oscillators 

We thus proceed to show analytical results of the en- 
tanglement dynamics in a second model, for which it is 
possible to understand the behavior of the entanglement 
in the presence of an effect which has no classical ana- 
log, namely, the quantum interference. This is done for 
a model of two resonant oscillators coupled nonlinearly 
which in certain regimes describe a coupled Bose-Einstein 
condensates (BEC). The terms of the Hamiltonian 119|) 
are 



k = huj (4«fc + |) (fe = 1,2), 
H\2 = ( a[a,2 + a\a\ ) + f<?9 ( a\ai - 



a\a 2 



1 



2 (24) 



Such Hamiltonian has been applied to describe an ideal 
two-mode BEC in the equal scattering length situation 
[o"i| with internal Josephson coupling \jsj^. Recently, the 
entanglement dynamics between the modes for initially 
disentangled coherent states, was exactly calculated j33j . 
We reproduce here the analytical result for the RLE: 
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(25) 



FIG. 1: Mean long-time entanglement for initial coherent 
states centered in (a) a periodic orbit and in (b) a chaotic 
trajectory. The insets show the result for the short-time en- 
tanglement, which behaves like S = 0.05 t 2 , for all values of J. 
Upper curves correspond to larger values of J, which assumes 
the values 3.5, 6.5 and 10.5. In these calculations we used 
e — u = 1.0, G = G' = 0.35 and h = 1, in arbitrary units. 

At least two aspects are remarkable in this numerical 
result: (i) the short-time behavior, proportional to t 2 , is 
J-indepcndcnt and (ii) the plateau value increases with 
J. 

The aspect (i) has already been predicted in section 
llll but for a different class of Hamiltonians. The results 
shown in Fig^are, therefore, a numerical verification of 
(|16f> in a more general context, in which the Hamilto- 
nian operator H is not obtained from a classical function 
by means of a quantization based on bosonic coherent 



in which f3\^{i) = («i,2e tut cos Xt — ict2,ie lu>t sm\t) 



and a>i = (qi + ipA/v2h. Result (|25|) predicts that the 
system recovers totally its purity (disentangled states) at 



t 



7T 



(26) 



after passing by an intricate interference process (see |33j ] 
for details). 

Following the reasoning presented in Sec[J our argu- 
mentation will be supported by a classical entropy, which 
is defined within the classical theory of ensembles, ft de- 
scribes the classical correlation dynamics of initially sep- 
arable probabilities. This quantity has been shown to 
evolve very closely to its quantum counterpart, the re- 
duced linear entropy, in cases that are initially separable 
and in which intrinsic quantum effects are absent |l3j| . 
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Similar approaches may be found in |34| . We then intro- 
duce the classical reduced linear entropy (CRLE), defined 
in Rcf.[ni as follows 



Sd(t) = 1 - 



ti-k[P k (t)} 
tr fe [P fe (0)]' 



(27) 



in which the partial classical trace is given by tr^ = 
J dqkdpk, for k = 1,2. Pi 2 denotes the reduced proba- 
bility function obtained by tr2,i[P(i)], in which the joint 
probability P(t) is a function satisfying the Liouvillc 
equation 



dP 



= {H,P} 



(28) 



TL is the classical Hamiltonian function. The normaliza- 
tion in (|27|l is necessary to guarantee both an adequate 
dimensional unity and S c i(0) = 0. The connection be- 
tween quantum and classical worlds at t = is estab- 
lished by taking 



P(0) 



(^1^2 1 po \v2v1) 
2nh 



(29) 



where po is the pure density operator at t = 0. The 
normalization is chosen in such a way that trP = 1. The 
initial distribution l|29|) is the unique source of h in the 
CRLE dynamics. Thus, the semiclassical regime in the 
classical Liouvillian formalism is associated to strongly 
localized initial distributions. 

The numerical results for both quantum and classical 
reduced linear entropies are shown in FigEI The depen- 
dence of CRLE in time reflects the fact that classical 
probabilities also become correlated, i.e., initially inde- 
pendent probability distributions are transformed into 
conditional ones as the dynamical evolution takes place. 
In some cases, similarities between these quantum and 
classical quantities indicate that entanglement has in- 
deed a strong statistical character ^3]. Here, however, 
remarkable differences emphasize the role played by an 
intrinsic quantum effect, as is the interference. In fact, it 
has been shown that interference is the major responsi- 
ble for such flagrant differences, since there is no analog 
effect in the classical formalism |33| . 

Concerning the semiclassical limit, we see that quan- 
tum recurrences, and consequently, quantum interfer- 
ences, are postponed by the dynamics as H is made small 
as compared with the classical action \fli\ 2 + |/32| 2 . This 
is indeed predicted analytically by (|2l)|> . Accordingly, 
the resemblance between quantum and classical entropies 
tends to increase in semiclassical regimes. In this sense, 
these results allows us to consider the CRLE as the very 
semiclassical limit of the entanglement. This is an indi- 
cation that in a semiclassical regime, the occurrence of 
the entanglement is due solely to spreading effects of the 
wave function, i.e., statistical effects present also in the 
classical theory of ensembles. 

The entanglement dynamics of the system l|24[) 
presents the same qualitative general features of our 
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FIG. 2: (a) Entanglement between two Bose-Einstein con- 
densates and (b) the classical reduced linear entropy for the 
classical counterpart as a function of time. Higher plateaus 
were attained for smaller values of h. The parameter values 
used were u) = 1, g = 0.1, A = 0.2, qi = Pi — 1 and h assumed 
the values 0.1, 0.5 and 1.0 (arbitrary units). 



precedent example, namely: short-time behavior inde- 
pendent of H (see [35( for the analytical demonstration) 
and higher plateaus for more classical regimes. Actually, 
this scenario has been shown to be quite general for a 
wide class of nonlinear systems [3(|. 

These results lead us to separate the dynamics of the 
entanglement in two main time scales. The first one, 
namely the short-time scale, may be regarded as a classi- 
cal scale in which entanglement is determined essentially 
by the vicinity of the initial condition in classical phase 
space. Accordingly, it has been shown that this scale is 
connected with the Ehrenfest time for the system l|24|) 
|33l |. Then, this is a situation in which, although the 
subsystems are entangled, the quantum state still allows 
the correspondence between quantum expectation values 
and classical trajectories. 

The second scale, associated to long times (for the 
which the Ehrenfest theorem does not apply), contains all 
the allowed quantum effects, like interferences, and also 
indicates the fact that the entropy is extensive with the 
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number of pure states accessible to the dynamics. These 
remarkably different time scales and the extensivity of 
the entropy have already been mentioned in slightly dif- 
ferent contexts in |3j| and j3?| respectively. Also, other 
time scales have been established in [35| . 

All results shown here point out to the same impor- 
tant conclusion: in closed pure bipartite systems, entan- 
glement exists even in arbitrarily scmiclassical regimes. 
On the other hand, in a strictly formal classical limit 
(H = 0), entanglement is indeed expected to disappear 
completely, since classical points have no statistics asso- 
ciated. In fact, this is a crucial assumption for obtaining 
the Newtonian trajectories from the quantum formalism 
. But this is a peculiar mathematical limit which will 
never be accessible in physical world. In this sense, our 
results also indicate that the formal classical limit of en- 
tanglement is rather singular for closed systems. 

The analysis presented above concerns a theoretical 
framework for the classical limit of quantum mechanics 
in closed pure systems. It is a mathematical limit which 
ensures that quantum mechanics is a universal theory 
which recovers the classical results in a particular regime 
of parameters. In closed systems, the limit is asymptotic 
for expectation values but seems to be rather singular 
for the entanglement. However, this conceptual scenario 
does not match the real world of quantum open systems. 
In this sense, our results must be regarded as the start- 
ing point for a more fundamental analysis that takes into 
account the effects of the environment, to which every 
real system is coupled. As it is well known, decohcr- 
ence tends to wash out quantum superpositions and the 
entanglement is expected to disappear under such situa- 
tion. However, a rigorous test of such assertions requires 
an adequate entanglement measure for mixed states of 
continuous variable systems, which is still a very contro- 
versial issue. 



IV. CONCLUDING REMARKS 

We have shown by two numerical examples and by an- 
alytical calculations that entanglement does not decrease 



as the semiclassical regime is attained asymptotically in 
closed pure systems. In fact, the results pointing out to 
an increasing amount of entanglement as the equilibrium 
is attained. 

This behavior may be understood by noticing that 
scmiclassical regime avoids quantum interferences, but it 
is not able to eliminate the spreading of the wave packet. 
Even extremely localized coherent states perceive the lo- 
cal spreading caused by the classical flux of trajectories 
associated. Furthermore, since we note that the entropy 
is extensive with the number of pure quantum states of 
the density operator, the increase in the long-time entan- 
glement in the semiclassical limit becomes just a natural 
fact. 

Recently, it was shown that entanglement can always 
arise in the interaction of an arbitrarily large system in 
any mixed state with a single qubit in a pure state |38j |. 
This result stresses the role played by the subsystem pu- 
rity as an enforcer of entanglement, even in a thermody- 
namic limit of high temperatures. Here, we are concerned 
with a different kind of classicality, namely, that one at- 
tained by h — > 0. Furthermore, we have studied strictly 
the case where the subsystems are in initially in pure 
states. The common point in these investigations is that 
the regarded systems do not interact with the environ- 
ment. In this sense, unitarity seems to be a key word in 
such apparent contradicting semiclassical limits. Then, 
there is a remaining question to be answered: Is indeed 
decoherence able to provide the "expected" semiclassical 
limit for the entanglement? 
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